Diagonal charge modulation in the insulating La 2 _ x Sr x Cu04 
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We show that, due to the Dzyaloshinskii-Moriya and the XY anisotropies, the disordered diagonal 
spin spiral in the insulating La2- :E Sr a ;Cu04 generates a diagonal charge density wave (CDW) with 
the wave vector twice of that of the spin spiral. The amplitude of the CDW depends on values of 
the anisotropies, the doping level, and on the density of states at the chemical potential. Based 
on available experimental data we estimate that for 4% doping the amplitude of CDW is about 
1/10 of the doping level. We believe that this mechanism explains the CDW observed recently in 
Zn-codoped detwinned La2- a: Sr I Cu04. 



The three-dimensional antiferromagnetic Neel order in 
I^-zSr^CuC^ (LSCO) disappears at doping x as 0.02 
and gives way to the so-called spin-glass phase which 
extends up to x as 0.055. In both the Neel and the 
spin-glass phase the system essentially behaves as an 
Anderson insulator and exhibits only hopping conduc- 
tivity [l], Superconductivity then sets in for dop- 
ing x > 0.055. The incommensurate magnetic order 
has been observed at low temperature in clastic and 
inelastic neutron scattering. According to experiments 
in the Neel phase, the incommensurability is almost 
doping-independent and directed along the orthorhom- 
bic b axis [3J. In the spin-glass phase, the shift is also 
directed along the b axis, but scales linearly with dop- 
ing 0,Sij. Finally in the underdoped superconducting 
region (0.055 < x < 0.12), the shift still scales linearly 
with doping, but it is directed along one of the crystal 
axes of the tetragonal lattice ■ 

In the present work we discuss only the insulating spin- 
glass phase, 0.02 < x < 0.055. The theory for the insu- 
lating phase has been suggested in Refs. @, @, The 
theory has the following essential components: 

1) Duc to strong antiferromagnetic correlations, the 
minima of dispersion of a mobile hole are at points 
(±7r/2, ±7t/2) of the Brillouin zone, so the system can, to 
some extent, be considered as a two valley semiconduc- 
tor. The hole does not have a usual spin, but it possesses 
a pseudospin that describes how the hole wave function 
is distributed between two magnetic sublattices. 

2) At low temperature, each hole is trapped in a hydrogen- 
like bound state near the corresponding Sr ion, the bind- 
ing energy is about 10 — 15mcV and the radius of the 
bound state is about 10 A. 

3) Duc to the orthorhombic distortion of LSCO, the ma- 
trix elements t' a and t' b describing the diagonal hopping of 
the hole are slightly different, and this makes the b valley, 
(— 7r/2, 7t/2), deeper than the a valley, (ir/2,ir/2). So all 
the hydrogen-like bound states are built with holes from 
the b- valley. In what follows, we refer to these bound 
states as impurities. 

4) Each impurity creates a spiral distortion of the spin 
background in the orthorhombic b-dircction. The distor- 



tion is observed in neutron scattering. So the state at 
0.02 < x < 0.055 is not a spin glass, it is a disordered 
spin spiral. 

5) At the point of overlapping of bound states ("percola- 
tion" point) the direction of the spiral must rotate from 
the diagonal to parallel because the Pauli principle. Si- 
multaneously the superconducting pairing is getting pos- 
sible. Hence we conclude that x = 0.055 is the percola- 
tion point. 

Intrinsically this picture does not contain any charge 
ordering and this is qualitatively different from the stripe 
scenario |11| . Charge modulation in the spin spiral pic- 
ture is certainly possible, but this can only be a secondary 
effect that is due to the spin-orbit interaction. Stimulated 
by the recent discovery [l2T | of the charge modulation 
in La1.g5Sro.05Cuo.95Zno.05O4 we suggest in the present 
work a specific mechanism for CDW that is due to the 
Dzyaloshinskii-Moriya (DM) and the XY anisotropies. 

The spiral pitch. Calculation of the spiral pitch 
can be performed within the mean-field approximation. 
Mobile holes are trapped by Sr ions in hydrogen-like 
bound states ("impurities") . The ground state of 
the "hydrogen atom" is four-fold degenerate: (two-fold 
pseudospin) x (two- fold valley). The orthorhombic distor- 
tion lifts the valleydegeneracy, so all the impurities reside 
in the b- valley [l~3| . Impurity pseudospin interacts with 
the spiral distortion of the spin fabric Q . The interaction 
energy is \p2~gQb where Q is the wave vector of the spiral, 
Qb is the component along the orthorhombic b-direction, 
and the coupling constant is approximately equal to the 
antiferromagnetic exchange, g as J as lAQmeV. We set 
the tetragonal lattice spacing equal to unity so the wave 
vector Q is dimensionless. The pseudospin degeneracy 
is lifted as soon as the spiral is established, all pseu- 
dospins are aligned and the corresponding energy gain 
per unit area is —xy/2gQ. Here x is concentration of im- 
purities that is practically equal to doping. The elastic 
energy of the spin fabric deformation is p s Q 2 /2 . Here 
p s as 0.18 J is the spin stiffness. Thus the total energy is 
p s Q 2 /2 — x\f2gQb. Minimization with respect to Q gives 



Q = Qb = x . 



(1) 
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To fit the experimental data 0, 0, Q we need g = 0.7 J 
that agrees with the t — J model estimate, g w J. The 
presented mean-field picture does not address the stabil- 
ity of the state, broadening of the line due to disorder, 
topological defects, etc. These issues have been studied 
in Refs. @, [nj. The stability depends on the localiza- 
tion length (size of the impurity) that does not appear 
in the mean-field picture. However, as soon as we know 
that the disordered spiral state is stable, then the above 
description is correct. We would like to stress that the 
spiral picture does not necessary assume a static spiral. 
The spiral can be dynamic. In particular in a pure 2D 
system the spiral is dynamic at any nonzero temperature. 
In LSCO, due to anisotropies and a weak 3D coupling, 
the spiral becomes dynamic at a small finite tempera- 
ture, T ~ 20K. However, the absence of the static spiral 
at T > 2QK does not mean that the spiral is not there, 
it just becomes dynamic. The important components in 
the above picture are l)the hole binding, 2)the a-valley 
depopulation, and 3)thc height of the spin-wave dome 
Ecross observed in neutron scattering. Both the bind- 
ing energy and the valley anisotropy energy are about 
10-15meV [1, [l3|. The value of E cross depends on dop- 
ing, and for x = 0.03 — 0.05 it is also about 15mev [14| . 
Therefore the spiral description is valid up to character- 
istic temperature Th ~ 150K. 

Density of states. In the case of uniform doping the 
diagonal spiral (unlike the parallel spiral) always has a 
tendency towards charge modulation [l5j, |l6| . In the case 
of the disordered state, the problem of charge instability 
was resolved in Refs. [US El, assuming that at zero tem- 
perature all the Sr-hole bound states are filled and hence 
there is no room for compression. So, implicitly the pic- 
ture of energy levels shown in FigJTjA was assumed: all 
the bound states are below chemical potential. However, 
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FIG. 1: The diagram of impurity energy levels. A: Chemical 
potential is outside of the impurity band. B: Chemical po- 
tential is inside the impurity band. Eb ~ 10 — 15meV is the 
binding energy, and 8E is the width of the impurity band. 

this picture would imply the activation behavior of dc 
conductivity a oc cxp{— (AE/T)} while it is well estab- 
lished that the conductivity follows the 2D version of the 



Mott variable range hopping (VRH) formula [l[ 

a ~ eM-iTo/T) 1 / 3 } . (2) 

This implies that the chemical potential is within the 
range of impurity energies, as it is shown in FigQj3. Hence 
some bound states are unoccupied and this gives room 
for CDW built on the bound states. Let us denote the 
concentration of unoccupied bound states by 6x. It is 
well established that the hole doping level is pretty close 
to concentration of Sr ions, therefore, Sx <C x. On the 
other hand 6x is not that small because it is sufficient for 
VRH. It is reasonable to assume that 
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(3) 



This is the maximum possible relative amplitude for 
charge density modulation, there arc no more quantum 
states within the impurity band to develop a larger am- 
plitude. 

The characteristic VRH temperature Tq in Eq. (J3J) de- 
pends on doping and sample quality and generally de- 
creases when doping increases (and thus conduction be- 
comes easier). At 4% doping the data of Ref. [l| are 
well fit with Tq « 500 K [17J. Analyzing the curves 
from Ref. we have found T (x = 0.02) - 8000^, 
T (x = 0.03) - 2000A-, T {x = 0.04) - 300^. The 
temperature Tq is related to the 2D density of states G, 



G = 



13.8 



(4) 



see Ref. [18j . Here I is the localization length (the "Bohr 
radius" of the bound state). We set the Boltzmann con- 
stant equal to unity. The localization length is about 2.5 
lattice spacings Note that the density of states de- 
termined in this way is valid up to T ~ E b ~ 100 — 150A" 
(ionization of bound states) in spite of the fact that 
the VRH resistivity formula is valid only up to T ~ 
20 — 30K. The density of states is related to the width 
of the impurity band SE shown in FigQ] G = x/5E. 
From here we find values of SE at different doping lev- 
els: SE(x = 0.02) - 6meV, SE(x = 0.03) - 2meV, 
SE(x = 0.04) - lmeV. 

Spin anisotropies and generation of the second har- 
monics of the spiral. It is very convenient to use the 
tr-model notation. The energy density in this notation 
reads [n| 



^ (Vn(r)) 2 - V2gJ2& ' X ( e b ' S(r - p< 



[D 2 n 2 a {v)+T c nl (r)] 



(5) 



2c 2 

Here n, n 2 = 1, is the staggered field that describes spins, 
£i is direction of pseudospin of i-th impurity, £ 2 = 1, ef, 
is the unit vector along the orthorhombic b-axis. As we 
have already mentioned, we neglect here the impurity 
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size. We mention once again that what we call "the im- 
purity" is the occupied bound state. The last two terms 
in ([5]) describe anisotropies induced by the spin-orbit in- 
teraction [III [2^ . The anisotropies "want" to direct n 
along the b-axis. The DM vector is D m 2.5meV and the 
XY anisotropy ^/T~ c f=a 5meV. The spin wave velocity is 
cw V2J. 

As we described above, the two first terms in ([5]) gen- 
erate the spiral 



(0, sin(Q -r + ip), cos(Q - r + ip)) 



(6) 



Here we assume that the spins are in the bc-planc. 
However, we will argue below that the actual plane of 
the spiral is not important, moreover, spins can slowly 
fluctuate without any static spiral. The disorder and 
topological defects give random phase ip that broadens 
the line (iq| . however the broadening is a separate is- 
sue and we disregard it here. The XY-term in ([5|) is 
tx n\ = [cos(Q • r + ip)f -> ± [1 + cos(2Q • r + 2(p)] -> 
— <psm(2Q ■ r). Here we have assumed that ip -C 1. 
Thus, the interaction does not vanish after integration 
over space only if 



ip = Asin(2Q ■ r) 



(7) 



This is the mechanism for generation of the second har- 
monics in the spin pattern. Substitution of © in (JSJ) 
gives the following energy density 



(Q + ^'f - V2ga(Q + ip') 



2c 2 



<psin(2Q ■ r) , (8) 



where cp' = (e& • V)<^ and <r(r) = J2i^( r ~ r i) * s the 
density of impurities per unit area of the plane. In the 
case of distribution of states shown in Fig. [TJ^ the density 
of impurities is constant because there is no room for a 
density modulation, a = x. In this case minimization of 
([8]) with respect to the amplitude of the second harmonics 
of the spin spiral gives 



A = 



r c 



8c 2 Q 2 



(9) 



The value of A is small: A ~ 10~ 3 — 10 -2 , so it is hardly 
possible to observe it directly in neutron scattering. 

Charge density wave. We know that the correct picture 
of states is shown in Fig. [T]B and in this case the CDW 
is possible. According to Eq. ((5]) the energy of a single 
impurity is shifted due to the second harmonics of tp as 



Se = - V2gtp' = -2V2gQAcos(2Q ■ r) 



(10) 



Hence, the variation of density of impurities is 
5a = -GSe 
B = 2GV2gQA 



Bcos(2Q ■ r) , 
13.8 T cPs 



0.3 r c 



8l 2 x J 2 T l 2 x JT 



(11) 



Here G is the density of states given by Eq. (|4|). Sub- 
stituting numerical values of parameters in (jlip we find 
that at x = 0.04 the amplitude of charge modulation is 



B ~ 10" 



(12) 



This value is in units of elementary charge per unit cell 
of the square lattice. Thus, the modulation is about 20% 
of the doping level. This dramatic enhancement of CDW 
is due to the very high density of states (|4]). In doing 
the estimate we assumed the static spin spiral in the bc- 
plane, see Eq. ([6|). However, since the effect comes from 
the phase (p, the static order is not essential. In the case 
of fully dynamic spiral we have to replace in (jlip r c — > 
i (r c + £> 2 ) . Since T c ~ D 2 , this does not influence the 
estimate (fT2|) . According to (fT2|) the CDW is so strong 
that it is quite possible that the real limitation on the 
amplitude B comes from the available Hilbcrt space limit 
given by Eq. ([3]). Anyway, both estimates (fT2")l and © 
give CDW amplitude of about 10% of the doping level. 

Interestingly, a naive selfconsistent treatment of the 
mean-field Eq. ((8|) gives an instability with respect to 
unlimited increasing of the CDW amplitude. However 
we know from the analysis 0, [l(| that the stability issue 
cannot be resolved within the mean-field approximation. 
This is why here we rely on perturbation theory. In any 
case, the amplitude is bound from the top by the condi- 
tion j3|) and since the perturbation theory result (fT2|) is 
of the same value as the upper bound we believe that this 
is a reliable estimate of the effect. The effect depends on 
temperature due to depopulation of the b- valley as well as 
due to ionization of bound states. The expected depen- 
dence is roughly B oc tanl^XK / '2T) . The characteristic 
temperature is T/, ~ 150-fT. 

Coulomb interaction, phonons, and large correlation 
lengths. The mechanism considered above explains the 
CDW amplitude, but it does not explain large corre- 
lations lengths observed in [12]. One needs an addi- 
tional weak interaction to coordinate the phase of the 
CDW. There are two candidates for the "coordination 
interaction", l)Coulomb interaction, 2)Interaction with 
phonons (lattice deformation). Let us first consider 
the Coulomb interaction. The dielectric constant k is 

For the direction along the c- 
30 - 70. The 



strongly anisotropic [21 1. 
axis it is equal to the ionic value k c 
in-plane value is much larger because polarizabilities of 
impurities contribute to the screening of the electric 
field [2l|. This contribution is proportional to doping, 
and extrapolating from data 2l| we get that at x = 0.04 
the value is K a b ~ 2000. Therefore, for estimates we will 
use the effective isotropic dielectric constant that is av- 
erage between k c and K Q b, K ~ 1000. Here we have in 
mind the zero temperature value. 

If the Coulomb interaction is important then it must 
establish the CDW antiphase between the CuC>2 layers. 
On the other hand we know that the in-plane modula- 
tion is very slow, 2tt/(2Q) <C d c , where d c = 13.15A 
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is the separation between the planes. A straightfor- 
ward electrostatic calculation shows that the Coulomb 
energy per unit area of the plane in this geometry is 
E c ~ j^3.2d c (e5a) 2 , where e5a is the charge density 
per unit area, e is the elementary charge and n is the di- 
electric constant (we use CGS units). The density wave 
is 5a = B cos(2Q • r). Therefore, with account of the 
electrostatic energy the impurity energy shift is changed 
from §M to 



5e = -2V2gQAcos(2Q-r) + 3.2—d c Bcos(2Q-r) . (13) 



The density variation is 5a = —G5e. Hence, we find from 
(|13|) that the density wave amplitude is 



B = 



2GV2gQA 
l + 3.2(e 2 /K)d c G 



(14) 



The difference from Eq. (jlip is in the Coulomb factor 
F c = [1 + 3.2(e 2 /K)4G] _1 . The value of this factor at 
x = 0.04 is Fc ~ 0.5. Thus, the Coulomb interaction 
does not qualitatively influence the estimate (fT5|) . It can 
change the estimate by at most a factor ~ 2. This con- 
clusion is supported by the following experimental ob- 
servation: The CDW amplitude is not very sensitive to 
temperature up to T ~ 100 - 150AT On the other 

hand, the screening must be very sensitive to tempera- 
ture because at T = the system is an Anderson insula- 
tor while at T > 50 — 70 A' it behaves like a conductor Q . 
Thus, the CDW amplitude is not sensitive to the change 
of the screening regime and hence the Coulomb interac- 
tion is not important. The most powerful confirmation 
of this point comes directly from experiment [T3] ■ In the 
observed CDW the layers are in-phase and this implies 
that the Coulomb interaction is negligible. 

Thus, as it has been pointed out in Ref 0, we are 
left with phonons to coordinate the CDW phase . Most 
likely this is also related to the DM interaction. The DM 
vector D is proportional to the oxygen octachedra tilt- 
ing angle, so we can write D — > D + 5D, where 5D is 
coupled to the soft phonon responsible for variation of 
the tilting angle. The £> 2 n^-term in ([5]) generates cou- 
pling to the spiral D 2 n 2 a —>(£) + <5D) 2 [sin(Q • r + ip)} 2 -> 
DSD cos(2Q • r + 2ip). This generates the lattice defor- 
mation at the second harmonics and this is the "coor- 
dination" interaction. Clearly, this mechanism gives the 
same phase of the CDW for nearest Cu02 layers. 

In conclusion, the disordered diagonal spin spiral in 
the insulating phase of LSCO generates a charge density 
wave with the wavelength half that of the spin spiral. The 
charge modulation is due to relativistic Dzyaloshinskii- 
Moriya and the XY anisotropics. At x = 0.04 we es- 
timate the amplitude of the modulation at ~ 10% of 
the doping level. The effect survives up to characteristic 
temperature 7\ ~ 150A". We believe that this theory 
explains the CDW observed in Ref. [l2j ■ 
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